Entanglement in the XX spin chain with an energy current 
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We consider the ground state of the XX chain that is constrained to carry a current of energy. 
The von Neumann entropy of a block of L neighboring spins, describing entanglement of the block 
with the rest of the chain, is computed. Recent calculations have revealed that the entropy in the 
XX model diverges logarithmically with the size of the subsystem. We show that the presence of 
the energy current increases the prefactor of the logarithmic growth. This result indicates that the 
emergence of the energy current gives rise to an increase of entanglement. 
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I. INTRODUCTION 

Recently entanglement properties of various quantum 
systems have been the focus of numerous studies. En- 
tanglement plays an essential role in several many-body 
quantum phenomena, such as superconductivity 1] and 
quantum phase transitions It is also regarded as an 
important resource in quantum computation and infor- 
mation processing Q. Quantum spin chains offer an 
excellent theoretical framevifork for investigating entan- 
glement properties, since several simple models can be 
solved analytically, and there also exist efficient numer- 
ical techniques. This motivated us to work with spin 
chains in order to investigate the effect of energy current 
on entanglement. 

There are two widely used method of characterising 
entanglement in spin chains. The first of these describes 
the entanglement between two spins in the chain with the 
quantity called concurrence 0, |3- The other one mea- 
sures entanglement of a block of spins with the rest of the 
chain with the von Neumann entropy, when the chain is 
in its ground state @, i7, 8j j3 , J.QJ . This latter method is 
especially useful when one tries to understand the role of 
entanglement in quantum phase transitions. These tran- 
sitions manifest themselves in the appearance of gapless 
excitations, and are accompanied by a qualitative change 
in the correlations. In view of the connection between en- 
tanglement and quantum correlations, the motivation to 
characterize a critical system in terms of entanglement 
properties naturally emerges. 

Vidal et al. calculated the von Neumann entropy for 
a wide range of one-dimensional spin models and found 
that for critical (gapless) ground states the entropy of a 
block of spins diverges logarithmically with the size of 
the block, while for noncritical chains it converges to a 
finite value. The prefactor of the logarithm was argued 
to be one-third of the central charge of the underlying 
conformal field theory. These results were supported by 



analytical calculations for the XX chain in and for 
more general Hamiltonians in j^] . 

Spin chains are simple enough models to investigate 
also the nonequilibrium effects on entanglement. One 
can find states that are characterised by the presence of 
currents of some phy sical quantities such as energy or 
magnetization [Tllll2j| . An important effect of these cur- 
rents is the rather drastic change in correlations. There- 
fore introduction of a current can be regarded as a quan- 
tum phase transition to a current-carrying phase. Con- 
sequently, it is interesting to find the entanglement prop- 
erties of these current-carrying states. 

In this paper we study an XX spin chain constrained 
to carry an energy current. We calculate the von Neu- 
mann entropy of a subsystem of L contiguous spins. The 
presence of the energy current maintains the logarith- 
mic asymptotics of the entropy; however, the prefactor 
of the logarithm is increased from 1/3 to 2/3, indicat- 
ing a higher level of entanglement in the current-carrying 
states. We also show that at a special value of the cur- 
rent, where the symmetry of the state is enhanced, the 
asymptotics of the entropy is the same as in the XX chain 
without current. In the vicinity of these transition points 
the entropy is shown to display a special type of finite-size 
scaling. 



II. XX CHAIN WITH ENERGY CURRENT 

The XX model is defined through the following Hamil- 
tonian: 
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where sf{a = x, y, z) are the Pauli spin matrices at sites 
Z = l,2,...,iVofa periodic chain and h is the magnetic 
field. Our aim is to constrain the spin chain to carry 
a prescribed amount of energy current; therefore we use 
the technique introduced in _12j . Since the local energy 
satisfies a continuity equation with the local energy cur- 
rent, one can calculate the operator of the total energy 



2 



current: 



N 



i+i 



i) 



1=1 



(2) 



h{s 



x„y 



In order to find the lowest-energy state among the states 
carrying a given current, one has to introduce a La- 
grange multipher A, and diagonahze the following modi- 
fied Hamiltonian: 
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The ground state of can be considered as a current- 
carrying steady state of at zero temperature. 

Since [H-^^ , J^] = 0, one can diagonalizc using 
the same methods which diagonalize H^-^ [l^, and the 
model can be transformed into a set of free fermions with 
the following spectrum: 



A/c — (— cos fc — — A sin k). 



(4) 



The ground state can be constructed by occupying all the 
modes with negative energy, and it remains the same as 
that of H^^ for A < 1. If the driving field A exceeds this 
critical value the energy current starts to flow, and the 
Fermi sea of the occupied modes splits into two parts. 
In order to illustrate the occupied regions it is useful to 
introduce the characteristic wavelengths kh = arcsin(/i) 
and k\ = arccos(A~^). The ground state can be analyzed 
as a function of h and the expectation value of the cur- 
rent density = {J^ /N), and the phase diagram shown 
on Fig^can be obtained. Three different phases can be 




FIG. 1: Phase diagram of the XX model with energy current 
in the h — plane, where h is the transverse field while 
is the density of the flux of energy. The occupied fermionic 
modes are represented by the black parts of the rectangles 

distinguished. In phase (g) and (3) only the magnetiza- 
tion energy part of the current is flowing, the current of 
interaction energy is zero, while the transverse magneti- 
zation, iW^, is nonzero. Entering phase the current of 



interaction energy starts to flow, while = through- 
out this region. On the line separating regions (T) and (f) 
(fc/i = kx) the symmetry of the ground state is enhanced, 
and it is characterised by a single Fermi sea. There are 
no states above the maximal current line, and in region 
@ the ground state is the same along the = const, x h 
lines, thus it can be represented by the h = 1 borderline, 
where the two Fermi seas merge. Details of the analysis 
of the phase space can be found in ■ 



III. ENTROPY OF A BLOCK OF SPINS 

We are interested in the ground-state entanglement be- 
tween a block of L contiguos s pin s and the rest of the 
chain. Following Bennett et al. |l4j we use von Neumann 
entropy as a measure of entanglement. It is defined as 

5l = -tr(/9L Inpi), (5) 

where the reduced density matrix pL = trjv-L|^g)(4'g| 
of the block is obtained from the ground state \^g) of the 
system by tracing out external degrees of freedom. 

In the calculation of the entropy we use a similar ap- 
proach that was succesfuUy applied in case of the XX 
model Q. The first step is to introduce the fermionic 
operators ci and c[ through the Jordan- Wigner transfor- 
mation. Note, that due to the nonsymmetric spectrum, 
we have to use fermionic operators instead of the Majo- 
rana operators that were used earlier in case of the XX 
model. The ground state in our case can be completely 
characterised by the expectation values of the two-point 
correlations (c]„c„) = Gmn] any other expectation value 
can be expressed through Wick's theorem. The matrix 
G reads 
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where the coefficients gi for an infinite chain [N 
are given by 
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where Ag is the spectrum defined in the previous section. 
Note that the integrand is just the characteristic function 
of the imoccupied fermionic modes. 

From the correlation matrix G one can extract the en- 
tropy Sl of the block as follows. First, by eliminating 
from G the rows and columns corresponding to spins that 
do not belong to the block, we obtain the correlation ma- 
trix Gl of the state pl- 
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In principle one can reconstruct the reduced density ma- 
trix pi^ using the matrix elements of G^. However, the 
entropy of the block Sl can be computed in a more di- 
rect way from the correlation matrix. Let U E S\J{L) 
denote a unitary matrix that brings Gl into a diagonal 
form. This transformation defines a set of L fermionic 
operators bm = X]n=i UnmCn which have a diagonal cor- 
relation matrix Gl = U~^GlU = diag(Ai, . . . , Al). The 
expectation values are thus 



0, {blb„ 



(9) 



that is, the above fermionic modes are uncorrelated. 
Therefore the reduced density matrix can be written as 
a product state 



PL = Pi 
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where /?„ denotes the mixed state of mode n. Hence the 
entropy pL is simply the sum of the entropy of each mode: 



asymptotically as ^ In L with the block size. As one starts 
to increase the value of the driving field A, the ground 
state (and the entropy as well) remains the same up to 
the critical field Ac = 1. Further increasing A one en- 
ters the current-carrying phase (2) (see Fignj, where the 
asymptotics of the entropy changes to |ln_L, and the 
same behaviour can be observed in phase . The only 
exception is the borderline of these phases (which is char- 
acterised by the condition = kx), where the entropy 
growth is again i In L. 

For values ft, > 1 of the magnetic field, all the spins are 
aligned in the ground state of the XX model, thus the 
entropy of a block vanishes. Switching on the current 
one observes a | In L entropy asymptotics in phase © . 

Summarizing the above results, one concludes that the 
introduction of an energy current may lead to a more 
rapid entropy growth, indicating a higher level of ground- 
state entanglement. The von Neumann entropy can be 
given for large block sizes (L ^ oo) as 
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A. Entropy asymptotics 

It follows from Eq. that in order to determine 

the entropy numerically, one only has to diagonalize the 
L X L matrix Gl, instead of diagonalizing the original 
2^ X 2^ reduced density matrix pL- This method reduces 
considerably the computational effort, and the entropy 
can be obtained for relatively large block sizes. Fig[21 
shows the results of the calculations. 
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FIG. 2: Entropy calculated from the reduced density matrix 
as a function of the block size L. The magnetic field is set 
\q h = 0.5; the curves correspond to difi'erent values of the 
current driving field A. The entropy grows as (2/3) In i in 
the current-carrying phases ® (A = 1.1) and @ (A = 1.3) 
except at the borderline (A = 2/-\/3), where the asymptotics 
is (1/3) InL, just as in the case of the XX model (A = 0). 

The ground state entropy of the XX model was first 
investigated in 0] and for /i < 1 it was found to grow 



S'l = ^ InL 
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where i? is the number of Fermi seas in the spectrum and 
5o is a function of the parameters h and A, independent 
of i. 

The above result was obtained analytically by Keat- 
ing and Mezzadri for general quadratic Hamiltoni- 
ans that have a correlation matrix of Toeplitz type with 
symmetric fermionic spectrum. In our case the correla- 
tion matrix is also of Toeplitz type, but the presence of 
the current breaks the left-right symmetry, resulting in 
a nonsymmetric spectrum. Nevertheless, the numerical 
results indicate that the above asymptotic form H12(l of 
the entropy seems to hold also in this more general case. 

The next to leading term Sq in the entropy is also 
given in a closed form in Q for symmetric spectra. Al- 
though the spectrum is not symmetric in our model, 
there are special cases when it can be transformed to 
a symmetric form. First we note that shifting the wave 
numbers by <y9 in the spectrum is equivalent to a unitary 
transformation V^GV of the correlation matrix, where 
V = diag(l, e*"^, e^*'^, . . . , e^^-^)*"^). Since this transfor- 
mation is diagonal, it leaves the eigenvalues of the re- 
duced density matrix Gl and thus the entropy invariant. 

Now, if hX = 1 (fc/i + k\ = ^) then the two intervals of 
the vacant fermionic modes (white parts of the rectangles 
on Fig^ have equal lengths, and a shift of the wave 
numbers by j symmetrizes the spectrum. In this case 
the constant term Sq can be expressed as follows: 
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where G = 1 -I- 7^ — 6/ In 2 is a constant defined through 
the Euler constant and I « 0.0221603 is a numerically 
evaluated integral expression Q. Thus one can see that 
the entropy can be written with a scaling variable as 
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FigEl shows the numericaUy calculated entropy with 
the logarithmic part substracted for different block sizes. 
The points perfectly fit the analytically calculated curve 
(|13|1 except near A^^ = which corresponds to the 
line separating the current-carrying phases (T) and (E) . 
Therefore, formula (fT!^ is applicable when £ 1. 
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FIG. 3: The next to leading order part of the entropy along 
the hX — 1 line for different block sizes as a function of A~^. 
The lines are analytically calculated from the form of the 
spectrum. The points perfectly fit the calculated curve except 
near the vicinity of = l/\/2 (see inset) and A^^ — 1, 
which corresponds to the boundary of the current-carrying 
phases. 

For nonsynimetric spectra the calculation of the next 
to leading order term in the entropy is mathematically 
more involved; hence we were only able to treat the prob- 
lem numerically. The results reveal that approaching the 
lines characterised hy kh = k\, kh — or kx — 0, Sq 
seems to diverge. This divergence is a consequence of 
the changing of the amplitude of the leading term: at 
the boundaries of the different phases the entropy grows 
as ^ In L. Thus approaching the boundaries, the | In L 
asymptotics has to be compensated with a negative log- 
arithmic divergence in 5*0. 



B. Finite-size scaling 

Obviously, for a fixed L, there is a finite neighbour- 
hood around the transition lines, where Ea. (|12|l cannot 
be used. In the case of a symmetric spectrum it was 
seen that it holds only when £ ^ 1. If we would like 
to characterise the behaviour of the entropy near these 
transition lines, we have to note that we can associate a 
diverging length scale, or alternatively a vanishing char- 
acteristic wave number to each of the lines. Similarly 



to finite-size scaling one writes the entropy near phase 
transition points as the sum of the "critical" entropy and 
a term depending only on the product of the block size 
and the characteristic wave number. For example near 
the high-symmetry transition line (kh = kx) it can be 
written as: 



SLikh,kx)^Sl + SiL\kh-kx\), 



(15) 



where S'£ is the value of the entropy on the transition 
line. The numerical calculations support the above type 
of scaling. Fig0| shows the numerical results for the scal- 
ing function near the high-symmetry line. Similarly, near 
the other transition lines (kx = or kh = 0) this type of 
scaling is valid, but with an other scaling function. 
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FIG. 4; Scaling function of the entropy near the kh = k\ 
transition line. The "critical' entropy is substracted from the 
total entropy and plotted against the scaling variable L{kh — 
k\). The curves are calculated for different lattice sizes (L = 
300 and 600) and for different values of the magnetic field 
{h — 0.3 and 0.7). The points all fit the same scaling function. 



IV. FINAL REMARKS 

We should note that analogously to the energy current, 
it is possible to introduce a current of magnetization. 
The resulting spectrum can be written in the same form 
as that of the XX chain, however with a shift in the wave 
numbers and a decreased effective magnetic field [l^ . 
Hence, the asymptotics of the entropy will be the same; 
only the L-independent constant increases. Thus, inter- 
estingly, the magnetization current has a much smaller 
effect on entanglement. 

In summary, we have shown on the example of the XX 
chain that the introduction of an energy current results 
a more entangled state. It would be worth considering 
current-carrying steady states in other spin models, to 
check whether the increase of entanglement is a general 
consequence of currents. 
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